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Abstract 

The Taylor coefficients of flavor diagonal and off-diagonal susceptibilities as well as baryon num- 
ber, isovector and electric charge susceptibilities are considered within a phenomenological quasi- 
particle model of the quark-gluon plasma and successfully compared with available lattice QCD 
data up to fourth-order for two degenerate quark flavors. These susceptibility coefficients represent 
sensible probes of baryon density effects in the equation of state. The baryon charge is carried, in 
our model, by quark-quasiparticle excitations for hard momenta. 
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I. INTRODUCTION 



The last few years witnessed two important milestones in the realm of relativistic heavy- 
ion collisions and related applied QCD: (i) Hints for a strongly coupled quark-gluon plasma 
have been deduced from experiments at RHIC [1-3], and (ii) lattice QCD calculations have 
been extended to non-zero net baryon density. While for many observables, at mid-rapidity, 
baryon-density effects are fairly small in heavy-ion collisions at top RHIC and future LHC 
energies, they become important for the ongoing low-energy runs at RHIC, previous CERN- 
SPS and future FAIR energies. Furthermore, the debated QCD critical point seems to be 
located at non-zero net baryon density according to investigations reported in [4]. There- 
fore, the exploration of this part of the phase diagram of strongly interacting matter gains 
increasing attention, both experimentally and theoretically. A necessary prerequisite in the 
search of the critical point is the understanding of thermodynamic bulk properties of QCD 
matter at non-zero net baryon density. 

First-principle lattice QCD calculations include all features of the complexity of QCD at 
finite temperature and net baryon density, supposed the numerical accuracy is appropriate. 
Indeed, signals of the QCD critical point have been found [5, 6], and the pseudo-critical 
curve not too far from the temperature axis in the T — \i plane is routineously determined 
today [5, 7] (T and /j, denote temperature and chemical potential, respectively). 

Basically, the partition function Z(T,/j,) or the grand canonical potential £l(T, //) or the 
pressure p(T, fi) contain much information on thermodynamic bulk properties of a medium. 
Susceptibilities are second-order derivatives of the pressure in the chemical potential di- 
rection. In so far, susceptibilities represent sensible quantities probing the active baryonic 
degrees of freedom in a medium. Even more, susceptibilities are related to fluctuations, 
which are debated to represent signatures of deconfinement effects [8, 9], thus being of 
utmost experimental relevance. 

Information on various susceptibilities have been accumulated from first-principle lattice 
QCD calculations [6, 10-17]. Keeping in mind possible limitations due to finite-size, nu- 
merical set-up and quark mass (m,) effects they, nevertheless, are a source of important 
information on baryon density effects in the hot quark-gluon medium. Below T c , where T c 
denotes the pseudo-critical temperature for deconfinement, the hadron resonance gas model 
(cf. [12, 18]) has been successfully compared with the lattice QCD data [12]. Above T c , the 
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situation is less settled. Certain baryonic bound states have been considered in [19] aiming 
at arriving at a physical picture of the strongly coupled quark-gluon plasma. Further devel- 
opments [20, 21] try to deduce also transport properties of deconfined strongly interacting 
matter. In addition, susceptibilities have been studied in phenomenological approaches such 
as the Nambu-Jona-Lasinio model [22] or Polyakov loop extensions thereof [23-25] as well as 
quasiparticle models [26, 27]. Furthermore, within the <fr functional approach to QCD [28], 
qualitative agreement with the lattice QCD data in [10] was found for T > 1.5 T c . All these 
approaches attempt to catch the relevant excitation modes. One should keep in mind that 
in the weak-coupling regime or in a medium with prominent quasiparticle excitations, the 
bulk properties are governed by excitations with hard momenta k ~ T, \i. Soft or ultra-hard 
modes are expected to influence p(T, fi) rather less. 

For the ultimate description of the very nature of the quark-gluon plasma one has to 
know correlations and spectral functions, propagators and related dispersion relations. Such 
information is still fairly scarce, but starts accumulating [29, 30]. Having at our disposal 
only the numerical data of thermodynamic state variables one must try to figure out which 
physical picture(s) is (are) compatible. Such an endeavor is clearly phenomenological. Be- 
sides the motivation of getting an interpretation of the available lattice QCD data, also the 
applicability of such phenomenological models represents important aspects, say for extrap- 
olations to larger net baryon density, or for interpolations between different regions of the 
QCD equation of state (EoS) [27], or for comparing different flavor numbers. 

The aim of the present paper is to confront in some detail the lattice QCD data from [6, 
12, 13] for two degenerate quark flavors with our quasiparticle model [31-33]. This model 
was used up to now for one independent chemical potential. With the goal of analyzing 
isovector, electric charge and flavor (off-) diagonal susceptibilities we are going to generalize 
the model towards two independent chemical potentials /i u and /i^ of up and down quarks, 
respectively. In fact, isovector and flavor (off-) diagonal susceptibilities represent much more 
sensible tests of a model than the baryon susceptibility alone. Furthermore, a detailed 
knowledge about the dependence of bulk thermodynamic quantities on, at least, two separate 
quark chemical potentials [i u and fi^ is necessary in order to discuss the impact of changes 
in various flavor sectors on the baryon density dependence of the EoS. Also, this becomes 
important when discussing properties of deconfined quark matter such as /^-stability and 
electric charge neutrality in hypothetical ultra-dense and hot proto-neutron stars. While for 
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one independent chemical potential, say fi u = fid, the model has been successfully compared 
with various sets of lattice QCD data in [33, 34], the straightforward generalization to a set 
of chemical potentials fl = {/i u , fid} is restricted by consistency requirements [35] given by 
Maxwell type relations and the stationarity condition of the thermodynamic potential. We 
show here, that the Taylor expansion coefficients of various susceptibilities are accessible up 
to a certain order in a consistent formulation, contrasting our model as an alternative to the 
picture developed, for instance, in [19]. 

Our paper is organized as follows. In section II, we extend the previous quasiparticle 
model [31-33] towards including two independent chemical potentials and discuss the con- 
sistency conditions for the resulting generalized system of flow equations. Section III is 
devoted to the numerical evaluation of various susceptibilities and the comparison with lat- 
tice QCD data. In addition, these results are used for discussing some properties of hot 
deconfined quark matter by means of a Taylor expansion of bulk thermodynamic quantities 
imposing, for instance, /^-equilibrium and electric charge neutrality. The summary and dis- 
cussion can be found in section IV. Appendix A summarizes the entropy density expression 
and its relation to the primary thermodynamic potential, while explicit representations of 
coefficients needed for determining the susceptibilities are listed in Appendices B, C and D. 



II. EXTENDING THE QUASIPARTICLE MODEL 

We consider the case of two degenerate quark flavors for which the lattice QCD 
data [6, 12, 13] are at our disposal. We choose the pressure p(T, fi U) fi d ) = ^lnZ(T, fi u , fi d ) 
as fundamental quantity in the following with quark flavor chemical potentials fi u> d or, equiv- 
alently, quark and isovector chemical potentials fi q j which are related via fi q = \{fi u + fid), 
fij = \{fi u — fid) or fi u = fi q + fij, fid = — Hi- (Note that fij was defined differently as 
either fii = 2(fi u — fid) in [11] or fij = \{fi u — Hd) in [6, 13].) fi q = and fii are associated 
with conserved quantum numbers (in strong interaction processes) of baryon number and 
isospin, respectively. The explicit model expression of p(T, fi u , fi d ) is relegated to Appendix 
A. The generalized quark number susceptibilities are defined by 



dfii u dfi ] d d 



(1) 

Ai«=Md=0 
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Because In Z(T, fi u , fi d ) is symmetric under CP transformations, derivatives for odd (j u + jd) 
vanish. Furthermore, in the flavor symmetric case m u = rrid = m, we find Xju,jd(T) — 

These generalized quark number susceptibilities represent a rich test ground. Besides the 
mentioned physical meaning of susceptibilities as measures for fluctuations, they addition- 
ally constitute the Taylor coefficients of the excess pressure Ap(T, fi u , /id) = p(T, fi u , fid) — 
p(T,fi u = 0, fid = 0), expanded simultaneously in powers of fi u and fid via 

A P (T,fi u ,fi d ) = J2x^(T)^$, (2) 

■ • Ju- 3d- 

thus containing information about baryon density effects in the EoS. (Similarly, one could 
consider the Taylor expansion of the pressure p(T,fi q + fii,fi q — fii) in terms of ji q and 
fij.) Various associated expansion coefficients, e. g. those of flavor diagonal and off-diagonal 
susceptibilities as well as relations among them, are discussed below. The excess pressure 
becomes increasingly important in the domain of larger values of fi U; d and lower temperatures. 
For instance, imposing /3-stability and electric charge neutrality of hot quark matter stars 
requires the knowledge about the dependence of bulk thermodynamic quantities on ji u and 
fid separately, at least. This underscores the importance of the susceptibilities Xj u ,jd(T) even 
if, at small T, a Taylor expansion in fi ud directions may not suffice. 
The net quark flavor number densities rii = dp/dfii with i — u,d read 

1 1 



di f°° 
rii = tt-tt / dkk 2 
Jo 



2tt 2 



e (ui-Hi)/T _|_ l e (ui+fn)/T _|_ i 



(3) 



where di = d Ut d = 2 N c refers to the spin and color degeneracies of the quarks. This implies 
that baryon charge- 1 carriers are quasiparticles with quark quantum numbers obeying the 
dispersion relations [36] 

u^k 2 + mt + Ui, U l = ^(r 2 + ^jG 2 (T,fi u ,fi d ). (4) 

For later purposes, we also exhibit the corresponding expression for gluons reading 

u; 2 g = k 2 + U g , u^^fr' + ^l + ^G'iT,^,^). (5) 

The quark mass parameters rrii might comply with the lattice calculational set-up, e.g. either 
rrii = with constant & to compare with [12], or constant m, to compare with [6, 13]. 
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The crucial point is that, besides the displayed explicit dependence of the self-energy parts 
III (where I is a label for u, d, g) on {/Xj} = /i Utd and T, there is also an implicit dependence 
via the effective coupling G 2 (T, (j, u , /j, d ). In the case of one independent chemical potential 
(J'u = = (i- e. /// = 0), the Maxwell relation dn q /dT = ds/d/i q , with entropy density s, 
together with the stationarity of the grand canonical potential, Sp/SUi = 0, leads to Peshier's 
flow equation [32] which determines G 2 (T,/i q ) for given initial condition G 2 (T,fi q = 0). In 
the case of two independent chemical potentials \i u ^ or, equivalently, /i q j, a system of three 
coupled equations is obtained from demanding stationarity and from the Maxwell relations 

ds dm 
ds dn q 

dfji q &r' [) 

dnj_ djh /ox 
dfi q d\ii ' 

The needed expression for the entropy density is listed in Appendix A and, with the defini- 
tions of jj, u> d and fi q j above, we note for isovector and quark number densities rij = n u — n d 
and n q = 3riB = n u + n d , respectively. The emerging system generalizes Peshier's flow equa- 
tion [32] towards two independent chemical potentials propagating G 2 (T,fi q = 0,/// = 0) 
into the thermodynamic parameter space, i. e. to non-zero ji q and ///. 

The structure of the generalized system of flow equations reads in the basis (fj, q , /ir) 

^£ +B ^ C2 ' (10) 

(A,-B,)— = (A 3 + B 3 )^, (11) 

with coefficients ^1,2,3, ^1,2,3, Ci,2 listed in Appendix B. From these coefficients it becomes 
evident how quark and gluon sectors are coupled. It was earlier argued [35] that the gener- 
alized system of flow equations in Eqs. (9-11) cannot be solved uniquely for arbitrary values 
of \x q and but only when assuming a side condition /i u = ii u (/i d ), i. e. when consider- 
ing merely one independent chemical potential. To see this, we reformulate Eqs. (9-11) in 
the basis (fj, u ,/i d ), making use of the analog of Eq. (11) in terms of ji u and fi d . Then, the 
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generalized system of flow equations is transformed into 



Ax 



A 



dG 2 

dG 2 



dG 2 



(12) 
(13) 



dpi u ' * dT 

These partial differential equations are uniquely solvable if the coefficients Aip and B\^ and 
Ci,2, as listed in Appendix C, are pairwise equal. Indeed, A\ = A2 and B\ = B2 hold in 
general. But, the furthermore needed equality C\ = C 2 is ensured for arbitrary but small 
values of fi Ut d, i- e. \i u ^ <C 7rT. (Actually, the equality of C\ and C 2 is given up to order 
0([i 2 ud ) in a Taylor series expansion in terms of /jl u and ji d \ the coefficients of third-order 
terms start to differ.) 

As the primary goal of this paper is the comparison of susceptibility Taylor series ex- 
pansion coefficients Eq. (1) (to be calculated as derivatives of p at fi Ut d = 0) with lattice 
QCD results up to fourth-order, the necessary conditions are fulfilled. The issue of potential 
limitations provided by the restriction to the small jj, u ^ region and one possible way of cir- 
cumventing them are discussed in Appendix D. In the needed leading order for evaluating 
the susceptibility coefficients of interest we note that from Eqs. (56) and (57) in Appendix C 

Ti 

i T 1 
-L-2 

X5 I\ 
X 4 X 2 

where Ik represent phase-space integrals listed in Appendix B, implying also u u = uj d in the 
mass symmetric case, m u = m d , and X 4 = X 5 . 



jJ<u = Vd: 



(14) 
(15) 



Furthermore, by exploiting Eqs. (9-11), one finds 



/iq=/i/=0 



dG 2 
din 



and 



8 2 G 2 



dp* 



flq=IH=0 



1 

AT 
+M 2 



2 £V + ! T G 2 (T) + l _ T 2^H 
w 3 V ; 3 dT 



(16) 



2e d T+ 2 -TG 2 (T)+ 1 -T 2dG2{T) 



while from Eqs. (12) and (13), we find |g 



dT 

2 ('n -I 5 ^G 2 (T) 
and 



Hu=fJ.d=0 



dG 2 



d 2 G 2 



1 

Jf 
-Is 



M 



' 2 e u T+ 2 -TG 2 (T) + l -T 2 ^l 



(17) 



1 G 2 (T)-I 4 ^- 2 G 2 (T) 



2tt 2 



with coefficients A/", A/1,2 listed in Appendix C, G 2 (T) = G 2 (T, /j, u = 0, fid = 0) and X 3i4i5 con- 
sidered at /i q — fj,i — or, equivalently, fJ> u — l^d — 0. Note that in the flavor symmetric case 
considered here, Eqs. (16) and (17) are related via d 2 G 2 /d/j, 2 \^ d=0 = \ d 2 G 2 /dfi 2 \^ ^_ Q . 
In addition, odd derivatives with respect to the chemical potentials such as or mixed 
derivatives such as a ^ or ^ ^ vanish at fi u — ^ — = [i q — [ij. The above stated ex- 
pressions and equalities are uniquely obtained from the generalized system of flow equations 
in the region of small ji u ^ and we can proceed by evaluating various susceptibilities. 



III. COMPARISON WITH LATTICE QCD DATA 
A. Taylor expansions in fi q /T at = 

In this section, we confront the above extended quasiparticle model (QPM) with lattice 
QCD data of various susceptibilities for Nf = 2 degenerate quark flavors. In [12], quark 
number and isovector susceptibilities as well as flavor diagonal and off-diagonal susceptibil- 
ities have been calculated on a lattice with temporal and spatial extensions N T = 4 and 
N a = 16 using improved actions and m u = ma = 0.4 T (i. e. £ u = £d = 0.4) as quark mass 
parameters. As a special case of the Taylor expansion in terms of \i q and jij, expansions in 
terms of fi q /T at /ii — were considered. 

The quark number susceptibility x q (T,H q )/T 2 = ^^0^ = 2c 2 + 12c 4 {^f + 

has been analyzed already in 

detail in [33], and an impressively good agreement of our model with the lattice QCD data 
from [12] has been found. The isovector susceptibility xi(T, fi q ) is only accessible with the 
present extension of our model; it obeys the expansion 

Xi(T,fi q ) d 2 (p(T,fi q + fi!,Hg- /i/)/T 4 ) j ! (H Q \ 2 , on / //V 4 



30c 6 (* ) 4 + OK) with c(T) = f^Z^T'^ 



T 2 d(fij/T) 2 
where the expansion coefficients read 

1 d k (T~ 4 p(T, ii q + fr, fx q - ///)) 



2cl + 12c{(^) + 30c^) +0(fi), (18) 



4(T) 



k\ d^j/T) 2 d(fi q /T) 



k-2 



(19) 



Due to the invariance of In Z under CP transformations, c{ vanish for odd k. 

For the first coefficient of interest we find within the extended QPM the explicit repre- 
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FIG. 1: (color online). Comparison of QPM results (solid curves) for cl with lattice QCD data [12] 
(circles for k = 2 and squares for k = 4) for iVj = 2 quark flavors. 



sent at ion 

a r°° u 2 p e ° 

*<r> = i?L dk T>(^rw> < 20 > 

with d = d u = da and e = uj u /T\^ u d=0 = u d /T\^ u d=0 which implies c\ = c 2 for all tempera- 
tures (cf. [33]) in the QPM. The next non-zero coefficient is 

= t4 / d^j-^—A e 2 - - 4e* + 1 (21) 



12tt 2 J T 3 (e e ° + 1) 4 

^ - 1) / 1 „ 2 ^ , T 2 «9 2 G 2 
^G (i ) + — 



e \ vr 2 6 <9/x 



where is given in Eq. (16). 

As in our previous studies [33], we choose for the effective coupling G 2 (T, \i u = 0, fi d = 0) 
entering Eqs. (4), (5) and (16) the parametrization 

, G 2 21oop (0, C = A^, T>T C , 
G\T) = { (22) 

The numerically evaluated QPM results for c 2 and C4 are exhibited in Fig. 1 (results for c 2j 4 
are exhibited in [33]) and compared with lattice QCD data [12], where we use as parameters 
entering G 2 (T) in Eq. (22) A = 5.95, T s = 0.75 T c and b = 421.5. The explicit value of T c is 
not important for the scaled quantities considered here. (Note that these parameters differ 
from the parametrization employed in [33] for describing q. This is due to the different 
quark dispersion relations used in [33] and here. Employing instead Eq. (4) as quark dis- 
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persion relation with the QPM parameters for G 2 (T) stated above, we find an equally good 
agreement of QPM results for q with the lattice QCD data [12] as reported in [33].) 

Similar to c 4 , the expansion coefficient c\ slightly underestimates the lattice QCD data 
[12] approaching its Stefan-Boltzmann limit 1 / (27r 2 ), while c\ agrees remarkably well with the 
data [12] for T > T c approaching its Stefan-Boltzmann limit Nf/2 asymptotically. Whereas 
Xq, due to the growing importance of the higher-order expansion coefficients with increasing 
chemical potential, exhibits a significant peak structure close to T c for large [i q jT indicating 
some critical behavior, xi does not point to such structures. This behavior is a consequence 
of the much less pronounced peak in c\ compared to C4. Similar findings were reported 
in [37], where a phenomenological sigma model was considered. Below T c , the agreement 
with lattice QCD data might be accidental, but one may consider Eq. (22) as convenient 
parametrization also for this region (see also discussion in section IV). 

Correlations between fluctuations in different flavor components can be discussed by 
considering flavor diagonal and off-diagonal susceptibilities. They read 

^ = \ + ¥) = 2 ° T + ucT ( t y + 3 ° cr ( t y + °^ (23) 

for the flavor diagonal susceptibility and 

^=1$-%)= ^ + (y) 2 + 3 ° c ^ ( t y + °^ (24) 

for the flavor off-diagonal susceptibility, where the individual expansion coefficients are de- 
fined by cl u = (c k + c[)/A and cf = (c k - c[)/A. 

The expansion coefficients c^ u and tff for k = 2,4 are exhibited in Fig. 2 and compared 
with lattice QCD data [12]. The diagonal expansion coefficients show a similar pattern 
as C2,4 and ^ approaching their Stefan-Boltzmann limits Nf/4 asymptotically in the case 
of c$ u and I /(An 2 ) for T > 2T C in the case of C4". The pronounced peak structure of the 
off-diagonal expansion coefficient c^ d is well reproduced, while in our extended QPM cVf 
is zero for all temperatures, in contrast to the data which are numerically small and differ 
noticeably from zero only in the region T < T c . This is simply a consequence of c\ = c<i- The 
pattern observed for the flavor off-diagonal susceptibility coefficients is discussed further in 
section III B. As the flavor off-diagonal susceptibility coefficients c^ d rapidly approach their 
Stefan-Boltzmann limit which is zero for all k, Xud vanishes for large T indicating that 
fluctuations in different flavor channels are uncorrelated at high temperatures. On the other 
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FIG. 2: (color online). Comparison of QPM results (solid curves) for the expansion coefficients c 
of the flavor diagonal susceptibility Xuu (left) and tff of the flavor off-diagonal susceptibility \ud 
(right) with lattice QCD data [12] (circles for k = 2, squares for k = 4). 

hand, Xud increases rapidly with increasing fi q in the vicinity of T c , indicating increasing 
correlations [8, 9] between fluctuations in different flavor channels in the transition region. 
This also explains the observed different behavior of Xq an d Xi'- While peak structures 
effectively add up in Xq they approximately cancel each other in xi- 

The behavior of the electric charge susceptibility xq is strongly related to Xq an d Xi v i& 
Xq — \ (xi + \Xq) ■ The corresponding Taylor expansion reads 

= 2c« + 12c? (^) 2 + 30c« (^) 4 + Oifi , (25) 

with expansion coefficients = | {d k + |cfe) = (|c^ u — |c^ d ). The coefficients and c? 
are exhibited in Fig. 3. From the definition of it becomes clear that contributions of 
pronounced structures appearing in flavor diagonal and off-diagonal susceptibilities do not 
completely cancel in the electric charge susceptibility. 



B. Taylor expansion in \i u and fid 

In this section, we confront the extended QPM with lattice QCD data [6, 13] of gen- 
eralized quark number susceptibilities Xj u ,j d as defined in Eq. (1). These simulations were 
also performed for Nf = 2 degenerate quark flavors on a lattice with temporal and spatial 
extensions iV T = 4 and N a = 16. However, the used quark mass parameter entering the 
quark dispersion relation reads now m u = = 0.1 T c , in agreement with the lattice perfor- 
mance [6, 13], which is temperature independent in contrast to the lattice set-up considered 
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FIG. 3: (color online). Comparison of QPM results (solid curves) for the electric charge suscepti- 
bility coefficients with lattice QCD data [12] (circles for k = 2, squares for k = 4). 

in section III A. As a result, some of the coefficients in the generalized system of flow equa- 
tions render which changes also the derivative expressions of the effective coupling. To be 
precise, the terms explicitly depending on £ u and ^, which enter Eqs. (16) and (17) and 
some coefficients in Appendix B and C, have to vanish for constant m u ^. 

Furthermore, non-improved actions have been employed in [6, 13], thus cut-off effects on 
the numerical results are sizeably increased compared to improved actions. In section III A, 
we assumed the lattice QCD data [12] to be rather close to the continuum limit as improved 
actions were used (cf. a discussion in [38]); thus no continuum correction factor was applied. 
(As discussed in [11], continuum limit corrections to the Taylor expansion coefficients Ck are 
expected to be similar (10-20%) to corrections for the pressure at zero chemical potential [39], 
even though, the corrections seem to increase for higher-order expansion coefficients, see [11, 
40].) Here, however, we have to rely on an estimate for the continuum extrapolation of 
the lattice QCD data from [6, 13]. By investigating different temporal lattice extensions 
N T at fixed large temperature in [41, 42], the continuum limit of some generalized quark 
number susceptibilities was estimated. Even though, in principle, correction factors could 
be different for different temperatures, we apply as scaling factors d^^ = 0.47 in the case of 
X2,o/T 2 [41, 42] and a larger correction = 0.32 in the case of X4,o [42] to the data [6, 13] 
for all T. 

Estimating the continuum limit is necessary for making possible a meaningful comparison 
between the expansion coefficients considered in section III A and the generalized quark 
number susceptibilities Xju,j d - m f ac ^ they are closely related [19], e. g. the expansion 
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coefficients of flavor diagonal and off-diagonal susceptibilities Xuu and Xud can be expressed 
in terms of Xj u ,j d via 

mu 1 X2,0 



uu 
C 4 



Aid 



2 T 2 ' 
l Xi,i 
2 T 2 ' 

2| (X4,o + 2X3,1 + X2, 2 ) , 
1 



24 



(2x3,1 + 2x2,2) 



Within the extended QPM, we find from Eq. (f) and by using Eqs. (28) and (29) 
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Xi,i(T) 
X4,o(T) 



d 



7r 2 , 



o, 

d 



IT 2 . 

71 ./o 
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3 2e 
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(26) 
(27) 
(28) 
(29) 

(30) 

(31) 
(32) 



(e 2e ° - 1) / 1 



-G 2 (T) + — 



T 2 <9 2 G 2 
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where 

8^ 



and |f 



M«,d=0 



are given in Eqs. (16) and (17). As both derivatives of 



the effective coupling entering these expressions are related with each other in the flavor 
symmetric case, we find X3,i — for all temperatures in the QPM, while %2,2 is non-zero. 
Furthermore, xi,i — as vanishes for all temperatures, while c% d is non-zero as %2,2 is 
non-zero. In particular xi,i, or c^, vanishes because flavor-mixing effects, which describe 
the dependence of one quark flavor sector on changes in another one, are inherent in the 
quasiparticle model only via the quasiparticle dispersion relations resulting in terms which 
vanish at fi u ^ = 0. Qualitatively, our findings, in particular the observed deviations in 
the flavor off-diagonal susceptibility coefficients, can be understood from perturbative QCD 
arguments. In a perturbative expansion of the thermodynamic potential different partonic 
sectors start to couple only at order 0(g z ) of the QCD running coupling g. However, these 
plasmon term contributions oc g 3 are not completely reproduced in a similar expansion of the 
quasiparticle model thermodynamic potential [43] . Similar findings, pointing to the necessity 
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FIG. 4: (color online). Comparison of QPM results (dashed curves for the parametrization em- 
ployed in section III A and solid curves for readjusted QPM parameters) for the generalized quark 
number susceptibilities X2,o/T 2 (left panel) and Xi,o (right panel) with the continuum extrapolated 
lattice QCD data [6, 13] (circles) and the lattice QCD data for 2c^ u from [12] (squares). 

of properly including flavor-mixing effects for affecting the flavor off-diagonal susceptibility, 
were reported in [44] within a Polyakov loop extended Nambu-Jona-Lasinio model approach. 

In Fig. 4, we exhibit the QPM results for X2,o/T 2 and X4,o and compare with the contin- 
uum extrapolated lattice QCD data from [6, 13] (circles). When using the QPM parameters 
found in section III A, the QPM results (dashed curve in the left panel of Fig. 4) underes- 
timate the lattice QCD data (circles) of X2,o/T 2 . For comparison, we also show the lattice 
QCD data [12] for 2c^ (squares), where the increasing deviations of the QPM results (dashed 
curve) from the data (squares) for increasing temperatures are due to the different quark 
mass parameters used here and in section III A. (Note that when applying continuum limit 
corrections of about 10% in the considered temperature range to the lattice QCD data [12] 
(squares) as stated above, both continuum extrapolated lattice QCD data sets [6, 13] (cir- 
cles) and [12] would be fairly well compatible apart from a narrow interval around T « T c 
such that one unique QPM parametrization would be sufficient.) To bridge the data for 
X2,o to X4,o by our model, we readjust, therefore, the QPM parameters entering G 2 {T) in 
Eq. (22) in order to perfectly describe the lattice QCD data [6, 13] (circles) of X2,o/T 2 by 
using A = 17, T s = 0.905 T c and b = 431. The corresponding QPM results for X2,o/T 2 and 
Xi,o are exhibited by solid curves in Fig. 4. Again, very good agreement for T > 0.9 T c is 
found. 
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C. Deconfined /3-stable and electrically neutral quark matter 



Now, we turn our attention to the discussion of some bulk properties of deconfined quark 
matter of Nf = 2 dynamical quark flavors by means of Taylor series expansions using the 
generalized quark number susceptibilities discussed in the previous section. Clearly, these 
considerations are limited by the range of validity of such an approach, say by conservatively 
guessing the quark flavor chemical potentials to be individually restricted by fi u ,d/T < 1. 
Note that we employ again m u ^ = 0.1 T c as quark mass parameter. 

Starting from the definition of the excess pressure Ap in Eq. (2), including only terms 
up to j u + jd = 4, the net baryon density n B , suppressing the explicit notation of the 
temperature dependence inherent in the generalized quark number susceptibilities, reads 

n B (Hu,Vd) = ^ | (X2,o + Xi,i ) (Vu + fJ-d) + (^ff + ^ff) (vl + Vd) ( 35 ) 
I 

+ 2 (X3,l + X2, 2 )(/4/^ + HuHl) 

n B (VB,Hl) = 77|2(x2,0 + Xl,l)^B + ^ ( 1X4,0 + ^3,1 + X2,2 ) 1^% (36) 



9 L 9 V3 ' 3' 

+ (X4,0 - X2,2)^BlA 



Thus, the net baryon density simultaneously depends on two independent chemical poten- 
tials, ji u and fid (or, equivalently, fi B and fij). This is similarly the case for a non- interacting 
gas of gluons and massless quarks with two independent quark flavor chemical potentials. 
Only in the special case of fii = 0, i. e. fi u = Hd = Hq, n B is a function of one chemical 
potential \i B alone ensuring constant net baryon density for constant baryo-chemical poten- 
tial. In general, however, a detailed knowledge about the dependence on different quark 
chemical potentials is required, when discussing baryon density effects on the EoS. This is 
illustrated in Fig. 5, where the scaled net baryon density is exhibited for constant ji B and 
constant temperatures. As by definition fi d = \\i B — fi u , one chemical potential in Eq. (35) 
can be replaced. We chose [i B jT = 1 such that fi u /T + fid/T = |, ensuring that these 
considerations stay within the range of validity of the employed Taylor expansion approach. 
The minimum at fi u /T = fid/T = | exhibits the value of n B /T 3 for one independent quark 
chemical potential. n B for \l b = T drops by 3.3% at T = 2T C and by 4.6% at T = 1.05 T c 
when changing /j, u /T from to |. Accordingly, one is tempted to consider the detailed 
knowledge about the individual fi u and jj, d dependencies as not so important. 
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FIG. 5: (color online). Scaled net baryon density tib/T 3 from Eq. (35) for constant Hb/T = 1 as 
a function of fi u /T for constant temperatures T/T c = 2, 1.5, 1.2, 1.1, 1.05 from top to bottom. 

However, there are physical situations, where the corresponding side conditions require 
the separate knowledge about the non-trivial \i u and fid dependencies of bulk thermodynamic 
quantities. First, we consider curves of constant which are given by the linear relation 
A*d = I/Ab — (see short-dashed curve in the left panel of Fig. 6 with /xs/T = 1). The 
individual net quark number densities read n u = X2,o^u + Xi^Hd + ^^H u ~^^'f J 'uf J 'd + ^-fJ' d + 
^f-HuH d and n d = X2,oHd + Xi,iA*u + + ^/^d + ir^u + ^vll^d- Since in the QPM 

Xi,i = X3,i = 0) lines of constant n u or are approximately given by lines of constant 
fi u or fid, i- e. simply vertical or horizontal lines in the left panel of Fig. 6 (not displayed). 
(Only at temperatures T « T c , where %2,2 is non-negligible, the simple pattern is deformed 
somewhat.) This situation is completely different when considering constant scaled net 
baryon densities as depicted by the solid curve in Fig. 6 (left panel) for ub/T 3 = 0.187 at 
T = 1.1 T c unravelling the non-trivial dependence of fid on fi u in contrast to constant fxs- 
In fact, here fis/T > 1 except for the case when fi u /T = or fidjT = 0. 

In heavy-ion collisions one often relates the quantum numbers of the entrance channel 
with the ones of the emerging fireball. Isospin-symmetric nuclear matter, for instance, is 
characterized by an electric charge per baryon ratio of 1:2. This translates into \n u — \n d = 
\tib which is fulfilled for fid = fi u , i- e. simply a diagonal line in the left panel of Fig. 6 (not 
displayed). Discussing, instead, gold on gold collisions, the electric charge per baryon ratio 
is approximately 0.4. The corresponding dependence fidi.Hu) for T = 1.1 T c is depicted by 
the dotted curve in the left panel of Fig. 6. Another important issue concerns electric charge 
neutrality in bulk matter. In pure Nf = 2 quark matter, electric charge neutrality would 
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FIG. 6: (color online). Left: Dependence Hdi/J-u) f° r various side conditions or physical situations. 
Hb = T is depicted by the short-dashed curve, whereas constant ne/T 3 = 0.187 holds along the 
solid curve where \ib > T. Electric charge neutrality is given along the long-dashed curve for 
pure Nf = 2 quark matter, while the dash-dotted curve includes additionally electrons, imposing 
/3-equilibrium. The dotted curve reflects the situation in Au + Au heavy-ion collisions. (The curves 
end where fx u /T + fXd/T > !■) Right: Scaled net number densities as functions of fi u /T demanding 
electric charge neutrality either for pure Nf = 2 quark matter (solid curves, n^/T 3 - top, n u /T 3 
- bottom) or for including electrons and requiring /3-equilibrium (dashed curves, n^/T 3 , n u /T 3 , 
n e /T 3 from top to bottom). For T = 1.1 T c . 

require |n u — ^nd = 0. The according dependence is depicted in Fig. 6 (left panel) by 

the long-dashed curve, again for T = 1.1 T c . More relevant for hypothetical very hot neutron 
star matter in a deconfined state is /3-equilibrium. Flavor changing weak currents give rise 
to the balance equation d <-> u + e + z/ e , i. e. in weak interaction equilibrium \x e = fid — f^u, as 
the produced neutrinos are supposed to leave the star and do not participate in the balance. 
The electron net density is approximated by n e = ^fi e T 2 + -^2^1, and electrically neutral 
bulk matter is determined by |n u — |n d — n e = 0. The corresponding dependence /id(/i u ) 
is depicted by the dash-dotted curve in Fig. 6 (left panel) for T = 1.1 T c . The d quark net 
number density decreases by requiring /5-equilibrium, demanding also a non-zero electron 
density for electrically neutral bulk matter (see Fig. 6 right panel); n u is not affected when 
including electrons and /3-equilibrium. This is in contrast to findings for the cold color-flavor 
locked phase of QCD [46] for Nf = 2 + 1 dynamical quarks, where no electrons are required. 
The discussion can easily be extended to the physically relevant case of two light (up and 
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down) and one heavier (strange) quarks, considering again two independent quark chemical 
potentials for the light quarks, fii = fi u = Hd, and for the strange quark, /i s . Recently, first- 
principle lattice QCD data for this case became available [38, 45]. A detailed comparison 
of the properly extended quasiparticle model with these lattice results and, in particular, 
a discussion of finite baryon density effects on the EoS relevant for the hydrodynamical 
description of the expansion stage of heavy-ion collisions demands further studies. 

IV. SUMMARY AND DISCUSSION 

The focus of the present paper is an analysis of isovector and various flavor (off-) diagonal 
susceptibilities for two-flavor QCD by comparing the extended quasiparticle model with 
lattice QCD data [6, 12, 13]. The model includes the same quark mass parameters m 8 as 
used in these lattice simulations. (Basically, one could accomplish also a chiral extrapolation. 
However, the effective coupling G 2 {T) may implicitly depend on these masses. This deserves 
separate investigations.) A crucial point to be kept in mind concerns finite size effects. 
The lattice QCD calculations [6, 12, 13] are performed on grids with finite temporal and 
spatial extension, while our phenomenological model is formulated in the thermodynamic 
continuum limit. The use of an improved action in [12] lets us hope that the finite size effects 
are sufficiently small to make a direct comparison meaningful. In contrast, the lattice QCD 
data of [6, 13] require severe continuum extrapolation factors. In so far, the comparison of 
our extended QPM with these data is less direct. 

Having these limitations in mind, we emphasize the good agreement of our model with the 
lattice QCD data for c 2j4 , c| 4 , c% u 4 and as well as for the related generalized quark number 
susceptibilities. We consider this successful comparison as encouraging. A conclusion is that 
quasiparticle excitations, with a mass gap also in the chiral limit, are able to explain those 
features of the strongly coupled quark-gluon medium which are encoded in the mentioned 
coefficients. In particular, baryon density effects are probed by these coefficients. The baryon 
charge is carried by quasi-quark excitations, in contrast to models [19], where di-quark and 
three-quark modes carry a substantial fraction of the baryon charge. Furthermore, in several 
physical situations, like relativistic heavy- ion collisions or in hot proto- (quark) neutron stars, 
the various mentioned coefficients are needed to implement the adequate side conditions. 

We have applied our model also for T < T c . Formally, the description of the lattice QCD 
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data below T c requires fairly large values of the effective coupling G 2 (T). (An alternative 
description could rely on strongly increasing correlations which are beyond the presently 
employed approach [47].) The corresponding excitations become very massive, ranging to 
hadronic mass scales. It turns out that a few massive excitations reproduce fairly well some 
of the lattice QCD data within the interval 0.8 T c — T c . This is numerically not too distinct 
from the hadron resonance gas model, where one may regroup several resonances into a 
few representative effective excitations. (Vice versa, we mention that the resonance gas 
model [11, 12, 48] coincides with lattice QCD data also slightly above T c ; for an even more 
extreme point of view we refer the interested reader to [49].) In this respect, it is conceivable 
that several models with fairly distinct assumptions may equally well reproduce the same 
lattice QCD data on thermodynamic bulk properties - examples are [19, 22-27, 50-52]. Only 
correlators and spectral properties of the excitations can unreveal their real nature in the 
strongly interacting system. 

On the other hand, the coefficient c^, and accordingly is poorly described. This 
may be a hint for missing modes or degrees of freedom in our model. Qualitatively, our 
findings can be understood since flavor-mixing effects, which are important for the correct 
description of the flavor off-diagonal susceptibility, are not explicitly inherent in our quasi- 
particle model, but only implicitly via the quasiparticle dispersion relations. Progressing 
lattice QCD calculations are welcome to resolve this issue and to get more confidence in 
the baryon number carrying modes (cf. discussions in [9, 18]). Also, the slight deviations 
between our model and the data very close to T c may signal a deficit of our quasiparticle 
picture. Nevertheless, considering our phenomenological model as useful parametrization of 
lattice QCD results, it may serve as QCD-based input for hydro dynamical calculations for 
the expansion dynamics of matter created in ultra- relativistic heavy-ion collisions, cf . [34] . 

Finally, we stress that the utilized Taylor expansion technique is sensitive to the region 
Hu,d — > 0. QCD critical point effects at larger values of fi u ^ may not be catched in such 
an approach. For a phenomenological procedure to supplement our model by critical point 
features see [53]. 

In summary, we extend our quasiparticle model towards two independent chemical po- 
tentials. This allows for the determination of various susceptibilities. We find an impressive 
agreement (with the exception of two numerically small flavor off-diagonal susceptibility 
coefficients) with lattice QCD data. Since a special set of susceptibilities also provides 
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the Taylor expansion coefficients of the baryon-driven excess pressure, we argue that our 
phenomenological quasiparticle model catches relevant modes for the equation of state at 
non-zero net baryon density. It may be used, therefore, for the future determination of 
higher-order Taylor expansion coefficients which become increasingly important at larger 
net baryon densities. 

Acknowledgements 

We gratefully acknowledge discussions with E. Laermann, F. Karsch, R. V. Gavai and 
S. Gupta. The work is supported by BMBF 06DR136, GSI-FE, EU I3HP. 

Appendix A 

The pressure p(T, fi u , fid) as primary thermodynamic potential of our model is constructed 
by assuming a quasiparticle picture via 



where B is to be determined as line integral from thermodynamic consistency conditions 
and the stationarity condition Sp/SUj = resulting in dpj/dUj = dB/dUj. The partial 
pressures pi of included excitations / referring to u quarks, d quarks and gluons (g) read 



where the dispersion relations ui = u>i(T, /i u , fid) are given in Eqs. (4) and (5), q is +1 
(— 1) for fermions (bosons), di refers to the spin (polarization) and color degeneracies of 
quasiquarks and quasigluons reading d u = dd = 2 N c and d g = N 2 — 1, and fi g = 0. In 
this way, we count left-handed transversal quasigluons as anti-particles of the right-handed 
ones. These structures emerge from the underlying two- loop QCD <fr functional [54, 55] 
by imposing formal manipulations such as neglecting finite width effects in the considered 
asymptotic HTL approximations of the one-loop self-energies, and neglecting (anti)plasmino 
and longitudinal gluon excitations as well as Landau damping [47]. While p is highly non- 
perturbative with respect to the effective coupling G 2 entering the self-energy expressions, 
it is this phenomenologically introduced coupling which enables the model to go beyond the 
^-derivable approximations in [55]. 




(37) 



l=u,d,g 



Pi = eidiT 
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The entropy density expression entering the generalized Peshier equations in Eqs. (6-8) 
reads s = Y.i= u ,d,g s i = d P/ dT with 



d 

-2d 



d 3 k 



(2tt) 3 
f d 3 k 
9 ' J^f 



ln(l + 



-(u>i-iXi)/T 



) + 



(ui - Hi)/T 



( e (ui-Hi)/T _|_ ]_) 



+ (fit 



In (1 - e-^' T ) - — 



UJg/T 



(39) 
(40) 



(e^/T- i)_ 

where i — u,d and /i u — /i q + fj,j, jid = /i q — fij. This additivity in the contributions s t of the 
various parton species is anchored in the underlying two-loop QCD <3> functional [47, 55, 56]. 
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The coefficients entering Eqs. (9 - if) read 



M = i 3 - 



B 1 



B 2 
C 2 



A, 
B 3 



+X 5 
-Xi 



2^ + ^W + ^) 



T 2 + -U q + ^f 



71 



+ x 2 



T 2 + ^(/i, -/x/) 5 



7T 



-Is—G [ii -l 4 - — [n q + n^G + J 5 -— [//g-/x/]G 
+1, U&T + 2 -TG 2 \ - J 2 f 2#T + jW 



2T 2 



2vr 2 



T 2 + —( N + f i I y 



IT 



3 



T 2 + ^ q -^f 



7T 
i 



T 2 + -( N + ^y 



71 



2*. 



T 2 + -u q -mf 



71 



1 12 12 

Z3 ^G 2 ii q -l 4 - — [n q + in] G 2 - J 5 - — [//, - // 7 ] G 2 

O 7T 



7T 



'3tt 2 

-i2 



T : i 2£ 2 T + ^TG 2 J + T 2 he d T + hc- 



I2 



1 



T 2 + -^ q + ^y 



71 



T 2 + -( N -^f 



(41) 



(42) 
(43) 

(44) 



71 



(45) 
(46) 

(47) 
(48) 
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where the phase-space integrals X^ are given by 

dn u d u f°° k 2 f e ("u+» q +»i)/T \ ( , q) 

1 dU u 27T 2 J 2U U T l v ( e K+^+Ml)/T + 1 )2 ( e ( Wu -M,-w)/T + 1)2 j ' ^ j 

dn d _ d d f°° k 2 ( e K+^-w)/ T e (<*d-i*+i*i)/T \ 

12 ~ dU~ d ~ 2^2 X 2^T V( e (- d +^-w)/r + 1)2 - ( e (- d -^+w)/T + 1)2 J ' <: ;)0 ) 

ds d f°° k 2 & )g l T 

X3= dn g = ^ 2 J dk 2T 2 (e^-l) 2 > (51) 

T _ ds u _ d u /■« A; 2 / K + /x g + ; u J )e^ + ^)/ r A 

J4 - sn: - _ 2^ y d *2^r* v — (e(-+M 9+w) /T + 1)2 — + uv - J A 52 ) 

Js = &n d = ~w I dk 2^F 2 { (e (^- W )/T + 1)2 + - J • ( 53 ) 

In Eqs. (14) and (15), /j g and /// in the phase-space integrals X fc are replaced by \i u = fi q + fii 
and n d = fx q - 
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The coefficients in Eqs. (12) and (13) read 



and 



A = A 2 
B 1 = B 2 

Ci = 



3 



4vr 2 
T 2 4- — 

9 



+ il 4 


T 2 1 


Ai 5 


T 2 + — 




9 
7T 2 


3 5 


9 



Zi ( 2#T + ^TG 2 ) - ls^G 2 fi u - T 4 ^ U G 2 , 



C! AwMr + 3 TG 



X, 



J_ r2(< (T 2 + ^/tt 2 ) 

2tT 2 ^X 2 (T 2 +/i 2 /7T 2 ) 



(54) 
(55) 

(56) 
(57) 



'37T 2 — X 2 (T 2 + /i 2 /^ 2 ) ' 

where /x g and /// in the phase-space integrals X& defined in Appendix B have to be substituted 
by /Ai = li q + /// and fi d = ^ q — Hi- The coefficients in Eq. (16) read 



M = 
M = 



2^ 



k 2 f e^ T n e 2 ^/ T 
T 2 u u V (e-«/ T + l) 2 (e^/ T + l) 3 



d n ' 



A: 2 



D 2uj d /T 



- 2- 



2tt 2 7 ~'"T 2 u d \{e^/ T + l) 2 (e^/ T + l) 3 
where X3 4 5 as well as uj Utd have to be taken at /i q — /ii — or, equivalently, /i u — /i d — 0. 



(58) 
(59) 
(60) 
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Appendix D 



Let us first briefly discuss an implication of the requirement [i u ^ <C ?rT needed for the 
consistency of Eqs. (12) and (13). Second-order susceptibility coefficients depend on G 2 
evaluated at /j, U:d = 0, while fourth-order coefficients depend on G 2 and d 2 G 2 /dfi 2 ld at 
Hu,d = 0. In general, n-th order derivatives of G 2 require up to and including (n — l)-st 
derivatives of C\ or C 2 - This implies that up to and including third-order the derivatives of the 
effective coupling can trustfully be taken. Therefore, second- and fourth-order susceptibility 
coefficients and related quantities are uniquely determined. However, 



and 

Mu=A 1 d=0 

higher orders cannot be evaluated uniquely. These derivatives enter, for instance, sixth- and 



higher-order susceptibility coefficients. 

The origin of this insanity is the special ansatz for the self-energy parts in the quasi- 
particle dispersion relations in Eqs. (4) and (5), while our primary thermodynamic potential 
in Eq. (37) together with (38) should allow for consistency in all orders of powers of jj, u ^- 
The reasoning for our ansatz in Eqs. (4) and (5) is the contact to one-loop expressions for 
the self-energies [57]. It has been shown, however, in [58], for one (imaginary) chemical 
potential, that one can discard the explicit /i 2 terms in the self-energies and obtains an 
equally suitable description of the lattice QCD results. In other words, the stationarity 
property of the thermodynamic potential p, involved in our quasiparticle model, causes a 
robustness against such modifications of the employed self-energy parametrizations. 

It happens that for the modified self-energies, IT = |T 2 G 2 (T, p, u , fi d ) and Tl g = 
|T 2 G 2 (T, /j, u , fi d ), the coefficients in Eqs. (12) and (13) render to 

T 2 

A 1 = A 2 = — (2J 3 + J 4 + J 5 ) , (61) 

B l = B 2 = ~X X , (62) 

Ci = C 2 = 2 -TG 2 Z x , (63) 

(for simplicity, we consider here the chiral limit or, as in section III B, temperature inde- 
pendent bare quark masses). I. e. the generalized system of flow equations in Eqs. (9-11) is 
uniquely solvable without restrictions, and G 2 and all its derivatives are trustfully obtained, 
implying also a consistent determination of the susceptibility coefficients of arbitrary order 
opening the avenue for future investigations. We have checked numerically that the result 
exhibited for the fourth-order coefficient in Fig. 4 is changed by less than 9% when changing 
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the self-energy expressions (generically a slight down shift of the curves occurs). The result 
for the second-order coefficient exhibited in Fig. 4 remains unchanged as it depends only on 
G 2 at ji u ^ d = which is not affected by the modification of the self-energy parametrizations. 
Similar statements are applicable for other related susceptibilities. Consequently, the results 
exhibited in Figs. 5 and 6 remain effectively unaltered. 
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